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Abstract 

We consider the Itzykson-Zuber-Eynard-Mehta two-matrix model and prove that the partition 
function is an isomonodromic tau function in a sense that generahzes Jimbo-Miwa-Ueno's [20]. In 
order to achieve the generahzation we need to define a notion of tau-function for isomonodromic 
systems where the ad-regularity of the leading coefficient is not a necessary requirement. 
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1 Introduction 

Random matrices models have been studied for years and have generated important results in 
many fields of both theoretical physics and mathematics. The two-matrix model 



was used to model 2D quantum gravity [14] and was investigated from a more mathematical 
point of view in [26, 17, 4, 5, 7, 12, 8]; the partition function of the model 



has important properties in the large A^-limit for the enumeration of discrete maps on surfaces 
[15] of arbitrary genus and it is also known to be a tau-function for the 2-Toda hierarchy. In the 
case of the Witten conjecture, proved by Kontsevich [24] with the use of matrix integrals not 
too dissimilar from the above one, the enumerative properties of the tau function imply some 
nonlinear (hierarchy of) PDEs (the KdV hierarchy for the mentioned example). On a similar 
level, one expects some hierarchy of PDEs for the case of the two-matrix model and possibly 
some Painleve property (namely the absence of movable essential singularities). The Painleve 
property is characteristic of tau-functions for isomonodromic families of ODEs that depend on 
parameters; hence a way of establishing such property for the partition function Z^^ is that of 
identifying it with an instance of isomonodromic tau function [20, 21]. 

This is precisely the purpose of this article; we capitalize on previous work that showed how 
to relate the matrix model to certain biorthogonal polynomials [26, 17] and how these appear 
in a natural fashion as the solution of certain isomonodromic family [9] . 

The paper extends to the case of the two matrix model the work contained in [9, 12, 10]; it 
uses, however, a different approach, closer to the recent [6]. 

In [9, 12, 10, 19] the partition function of the one-matrix model (and certain shifted Toplitz 
determinants) were identified as isomonodromic tau functions by using spectral residue formulcB 
in terms of the spectral curve of the differential equation. Such spectral curve has interesting 
properties inasmuch as -in the one-matrix case- the spectral invariants can be related to the 
expectation values of the matrix model. Recently the spectral curve of the two matrix model [8] 
has been written explicitly in terms of expectation values of the two-matrix model and hence 
one could use their result and follow a similar path for the proof as the one followed in [12]. 
Whichever one of the two approaches one chooses, a main obstacle is that the definition of 



d/i(Mi,M2) = e-T'-(^i(A^i)+^2(*^2)-Af,M2)d^^^d^2 



(1.1) 




(1.2) 
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isomonodromic tau function [20, 21] relies on a genericity assumption for the ODE which fails 
in the case at hand, thus requiring a generalization in the definition. 

According to this logic, one of the purposes of this paper is to extend the notion of tau- 
function introduced by Jimbo-Miwa-Ueno's [20], to the two-matrix Itzykson-Zuber model. This 
task is accomplished in a rather general framework in Sec. 3. 

We then show that the partition function has a very precise relationship with the tau-function 
so introduced, allowing us to (essentially) identify it as an isomonodromic tau function (Thm. 
3.4). 

2 A Riemann Hilbert formulation of the two- matrix model 

According to the seminal work [26, 17] and following the notations and definitions introduced in 
[4, 5], we consider paired sequences of monic polynomials {TTm{x),crm{y)}m=Q...oo {rn = degvTm = 
degcJm), that are biorthogonal in the sense that 

dxdy7r^(:E)cT„(y)e-^i(")-^2fe)+xj; ^ ^^^^^^ q (2.I) 

The functions Vi{x), V2{y) appearing here are referred to as potentials, terminology drawn from 
random matrix theory, in which such quantities play a fundamental role. 

Henceforth, the second potential V2{y) will be chosen as a polynomial of degree (i2 + 1 

d2+l 

V2{y)=Yl -y'' ^d,+i /O (2.2) 

For the purposes of most of the considerations to follow, the first potential Vi{x) may have very 
general analyticity properties as long as the manipulations make sense, but for definiteness and 
clarity we choose it to be polynomial as well. 

The symbol J stands for any linear combination of integrals of the form 

/ / dxdy := dx dy, Xij G C (2-3) 

J J X j k 



where the contours {Tk}k=i...d2 will be chosen as follows. In the y-plane, define ^2 + 1 "wedge 

2k 

d2- 



sectors" {Sk}k=o...d2 such that Sk is bounded by the pairs of rays: := {y\ argy = 9 + 



and := {y\ argy = 9 + ^^j^-|frp-}, where 9 := argt^^^j+i- Then f ^ is any smooth oriented 

contour within the sector Sk starting from 00 asymptotic to the ray rj. (or any ray within the 
sector that is at an angle < 2(2^2+1) ^^^^^ equivalent for purposes of integration), and 

returning to 00 asymptotically along r^-i (or at an angle < 2(2^2+1) ^^)' These will be referred 
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to as the "wedge contours". We also define a set of smooth oriented contours {f fc}fc=i,...d2' that 
have intersection matrix Tj n = 6jk with the Tfc's, such that starts from oo in sector Sq, 
asymptotic to the ray fo := {y\ arg(y) = 6 — and returns to 00 in sector Sk asymptotically 

along the ray ffc := {?/| arg(y) = 6 -\ — gj+f-- These will be called the "anti-wedge" contours. 
(See Fig. 1.) The choice of these contours is determined by the requirement that all moment 
integrals of the form 

/ y^e-^^^y^+'^ydy, [ /e^^fe)-^^^^^ k = l,...d2, jGN (2.4) 

be uniformly convergent in x G C. In the case when the other potential Vi{x) is also a polynomial, 
of degree di + 1, the contours {Tk}k=i,...di in the x-plane may be defined similarly. 
The "partition function" is defined here to be the multiple integral 

TV N 

^^ '■=m n dx.dy, A(x)A(y) J] ^-V,{xo)-V2{y,)+x,y, (2.5) 

^^■JJx^j^i i=i 

where A(X) and A(y) denote the usual Vandermonde determinants and the factor is chosen 
for convenience. 

Such multiple integral can also be represented as the following determinant 

Zn = det[/i,,]o<i,,<7V-i , := [ x^y^e-^i^^'^-^^^^'^+^Mxdy (2.6) 

The denomination of "partition function" comes from the fact [26, 17, 9] that when x coincides 
with M X M then Z]\f coincides (up to a normalization for the volume of the unitary group) with 
the following matrix integral 

J j dMidM2e-*'^(^i(*^i)+^2(*^2)-MiM2) (2.7) 

extended over the space of Hermitean matrices Mi, M2 of size N x N, namely the normalization 
factor for the measure dfi{Mi, M2) introduced in 1.1. 

2.1 Riemann— Hilbert characterization for the orthogonal polynomials 

A Riemann-Hilbert characterization of the biorthogonal polynomials is a crucial step towards 
implementing a steepest-descent analysis. In our context it is also crucial in order to tie the 
random matrix side to the theory of isomonodromic deformations. 

We first recall the approach given by Kuijlaars and McLaughin (referred to as KM in the 
rest of the article) in [25], suitably extended and adapted (in a rather trivial way) to the setting 
and notation of the present work. We quote -paraphrasing and with a minor generalization- 
their theorem, without proof. 
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Figure 1: Wedge and anti- wedge contours for V2{y) of degree D2 + I 



Theorem 2.1 (Kuijlaars and McLaughin asymptotic). The monic bi-orthogonal polynomial 
TTnix) is the (1,1) entry of the solution T{x) (if it exists) of the following Riemann-Hilhert 
problem for T{x) . 

1. The matrix T{x) is piecewise analytic in C\ 

2. the (non-tangential) boundary values ofr{x) satisfy the relations 



r(x)+ = r(x)_ 



1 Wj^l ... Wj^d2 

1 



d2 



^2 p 

k=i "'^fc 



(2.8) 



(2.9) 



3. as x ^ 00 we have the following asymptotic expansion 





r{x)^{ia + ^ + o 



1 



x" 



^-'"^-iM 














x-'^^Idd. 



■2-rNj 



(2.10) 
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where we have defined the integers rriN^r^ as follows 



N = mNd2 + rN, mN,r]\r£N, < r„ < (i2 



(2.11) 



It follows from [25] that the solution T]sr{x) has the following form 



Tn{x) := T{x) :-- 



C:{f{z)) :-- 



1 

2TTi 



PN-l{x) Co{pN~l) ■ 
PN-diix) Co{pN~d2) ■ 



Cd2-i(^Jv) 

Cd2-l{PN-l) 
Cd2-l{PN-d2) 



(2.12) 



(2.13) 



where the polynomials denoted above by pn-i, ■ ■ ■ ,PN-d2 s-re some polynomials of degree not 
exceeding — 1, whose detailed properties are largely irrelevant for our discussion; we refer to 
[25] for these details. 

By a left multiplication of this solution by a suitable constant matrix we can see that the 
matrix 

TTn Co{TTn) ■ ■ ■ Cd2-li'n'n) 

VTn-l Co(7r„_i) ... Q2-l(^n-l) 



N 



T^n-d2 Co(vr„_rf2) 



Cd2~l{'^n~d2) 



and Fjv are related as 



Tm{x) = Un^n{x) 



(2.14) 



(2.15) 



where Un is a constant matrix (depending on N and on the coefficients of the polynomials but 
not on x). As an immediate consequence, solves the same RHP as T with the exception of 
the normalization at infinity (2.10). 

The present RHP is not immediately suitable to make the connection to the theory of 
isomonodromic deformations as described in [20, 21]; we recall that this is the theory that 
describes the deformations of an ODE in the complex plane which leave the Stokes' matrices 
(i.e. the so-called extended monodromy data) invariant. The solution Tn (or Tn) does not solve 
any ODE as formulated, because the jumps on the contours are non constant. If -however- we 
can relate Ttv with some other RHP with constant jumps, then its solution can be immediately 
shown to satisfy a polynomial ODE, which allows us to use the machinery of [20, 21]. This is 
the purpose of the next section. 
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2.2 A RHP with constant jumps 



In [9] the biorthogonal polynomials were characterised in terms of an ODE or -which is the same- 
of a RHP with constant jumps. In order to connect the two formulations we will use some results 
contained in [13] and we start by defining some auxiliary quantities: for 1 < k < d2, define the 
d2 sequences of functions {'ipm\x)}meN as follows: 



and let 



1 
27ri 



ds 



dzdw '^'^^''^'' ^'^'^ ^^^""^ ~ c-^2{u,)+V2(s)+zu,-xs 



X — z 



w 



I<k<d2, 
(2.16) 

(2.17) 



In terms of these define, for > ^2, the sequence of {d2 + 1) x (d2 + 1) matrix valued 

... #)(x) 



functions $(x 

N 



N 



i;^^\x) 



The following theorem is easily established using the properties of the bilinear concomitant and 
it is a very special case of the setting of [3] (Cf. Appendix B for a self-contained re-derivation) 

Theorem 2.2 (Jump discontinuities in \I/). The limits when approaching the contours T,- 

N N 

from the left (+) and right(—) are related by the following jump discontinuity conditions 



(2.18) 



$+(x) = $„(3;)H(^) 
N N 



(2.19) 
(2.20) 



where 



H(^) := I - 27rieoK^ 



I -|- 2TTieQK 



(2.21) 
(2.22) 



The proof of this theorem is given in Appendix B. For later convenience we define also 

^ := UZ^^ (2.23) 

N N 





/1\ 




( ° \ 











eo := 




K := 






w 
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The relationship with the matrices Fat, Fjv introduced in the previous section is detailed in the 
following 

Theorem 2.3 (Factorization theorem). The following identities hold 



^(x) = Tn{x)V{x)W{x) , ^^^{x) = Tn{x)V{x)W{x) 
N N 



where 



V :-- 



-Vl{x) Q 



Vb, 

and Vq, Wo{x) are the d2 x d2 matrices with elements 



W{x) :-- 



1 

Wo{x) 



(2.24) 
(2.25) 



{Vo)jk 



(Wo(x)),fc 



V2 
t>3 



V3 



Vd2 + 1 



1Jd2 + l 



Vd2 Vd2 + 1 

Vj+k j + k <d2 + l 
if j + k> d2 + l, 



yi-ieV2(2^)-^s/^y^ l<j,k<d2 



(2.26) 

(2.27) 
(2.28) 



The proof is a direct verification by multiplication by matrices, noticing that the matrix 

Vq is nothing but the matrix representation of ^^iv) ^^2i^) ^ quadratic form in the bases 

y ^ 



l,y,y^,...,y'^^ ^ and l,s,s^, 



d2-l 



(more details are to be found on appendix A, based on 



[13, 3]) The RHP for vj/ can be read off from that of r^r and the fact that the jumps are constants. 

N 

For convenience we collect the information in the following 

Theorem 2.4. The matrix vj/ is the unique solution of the following RHP: 

N 

1. Constant Jumps: 



^+(x) = ^_(x)h(^) 

TV N 



(2.29) 
(2.30) 



2. Asymptotic at infinity: 



^(x) ~ Fat 

JV 



Np-Vi{x) 



' X e 








-miv — 1 



Id. 



rN 






x-""^Id. 



"^oix) 



'd2-rM/ 



(2.31) 
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where 

r7v = /d+^ + ... (2.32) 

X 

and where ^'0(2;) := V{x)W{x) will be referred to as the bare solution. Its asymptotic at 
infinity can be computed by steepest descent, but since it is N -independent, for the sake of 
brevity, we do not report on it (details are contained in [5, 13]). 

3. has constant jumps 

4. ^'^(x)^'^^ = Dn{x) where D]\f{x) is a polynomial in x 

5. duj^^ n{x)'^]^^ = Uk,n{x) is polynomial in x. 

6. d^j'^ n{x)'^J^ = Vj^Ar(x) is polynomial in x. 

7. det(^'Ar+i^']^^) = Cste 

The points (4,5,6,7) in the above theorem can be found in [8, 5] 

In the next section we shall define a proper notion of isomonodromic tau function: it should 
be pointed out that the definition of [20, 21] cannot be applied as such because -as showed in 
[5]- the ODE that the matrix ^ (or ^) solves, has a highly degenerate leading coefficient at the 

N N 

singularity at infinity. 

In the list, the crucial ingredients are the differential equations (in x or relatively to the 
parameters uk and vj). First, the fact that Dj^{x) is a polynomial comes from explicit com- 
putation (See [8] for example). The result concerning the determinant of Rn{x) can also be 
found in [8] where one has: det(^'Ar+i^']Y^) = det(aAr(x)) = Cste. The properties concerning 
the differential equations relatively to parameters can be found in [8] too. Under all these as- 
sumptions, we will show that the proof of Jimbo-Miwa-Ueno can be adapted and that we can 
define a suitable r- function in the same way Jimbo-Miwa-Ueno did it. 

3 Definition of the r-function 

In this section, we will place ourselves in a more general context than the one described above; 
we will show that under few assumptions one can define a good notion of tau-function. 

More generally we will denote with ta the isomonodromic parameters (in our case they are 
the uks and ths vfs) and a subscript a or 6 is understood as a derivation relatively to ta or t^. 
For a function / of the isomonodromic times we will denote by the usual symbol its differential 

df = Y,dtJdta = Y,fa<^ta (3.1) 

a a 
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Our setup falls in the following framework that it is useful to ascertain from the specifics of the 
case at hands. Suppose we are given a matrix 

^(x) - Y{x) E{x) , Y{x) := (l + ^ + ^ + ..}jx^ (3.2) 

where = E{x;t) is some explicit expression (the "bare" isomonodromic solution) and S is 
a matrix independent of the isomonodromic times. This implies that if we define the one-form- 
valued matrix H{x; t) by 

H{x; t) = dE;(x; t) E{x; t)~^ (3.3) 

then H{x) = J2'^a<ita (we suppress explicit mention of the t dependence henceforth) is some 
solution of the zero-curvature equations: 

daHb-dbHa=[na,Hb] (3.4) 

We will Eissume (which is the case in our setting) that all Ha are polynomials in x. We 
will also use that the dressed deformations fla given by = ^a* are polynomials. Moreover, 
according to the asymptotic they are given by: 

na = {YnaY-\oi- (3.5) 

In this very general (and generic) setting we can formulate the definition of a "tau function" as 
follows 

Definition 3.1. The tau- differential is the one-form defined by 

uj ■= ^Wadi" := ^restr {y-^Y'Ha) dt" (3.6) 

a 

The main point of the matter is that -without any further detail- we can now prove that the 
tau-differential is in fact closed and hence locally defines a function. 

Theorem 3.1. The tau-differential is a closed differential and locally defines a r-function as 

dlogr = a; (3.7) 

Proof. We need to prove the closure of the differential. We first recall the main relations 
between the bare and dressed deformations 

daY = QaY - YHa ; YHaY'^ = - 7^„ ; Ka := daYY'^ (3.8) 
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We note that -by construction- = {yHaY ^)poi is a polynomial while TZa = 0{x ^) irrespec- 
tively of the form of S. We compute the cross derivatives directly 



resTvJ^ - {^aY - YHa) Y'^Y'Tib + Y'^ {QaY - YHa)' Hb + Y-^Y'daHb^ 

lesTt (UaY-^Y'Tib + Y-^n'^YHb - Y^^Y'TiaHb - H'jib + Y-^Y'daH^^j 

polynomial 

resTv(^y-ir {[Hb^Ha] + daUb) + Y'^^'^YHb - ) 

resTY(^y-ir {[Hb^Ha] + daHb) - ^'a^^ (3.9) 



where, in the last step, we have used that YTLbY^^ = $1^ — IZb and that the contribution coming 
from ^Ib vanishes since it is a polynomial. Rewriting the same with a ^ b and subtracting we 
obtain 

dacob - dbUJa = TesTr(2Y-'Y'[nb, Ha] - ^TZb + ^'^,7^a + Y-^Y' {daHb - dbHa) ^ 

=0 by the ZCC 3.4 

= TesTi{Y-^Y'[nb,na\ - KTZb + n'^TZa + Y-'Y'CdaHb - dbHa + [^6, Wa] 



resTY (Y-^Y'iHb, Ha] - n'^TZb + l^x) 



(3.10) 



Note that, up to this point, we only used the zero curvature equations for the connection 
V = Y,{da - Ha)dt" and the fact that Ti a are polynomials in x. We thus need to prove that the 
last quantity in (3.10) vanishes: this follows from the following computation, which uses once 
more the fact that 7ia and are all polynomials. Indeed, we have resTr('H^'Hfe) = and hence 
(using (3.8)) 



= iesii{n'Jib) = resTV {{YnaY~^)' YHbY-^j - resTV (Y'TiaHbY-^) + resTV {HaY-^Y'Tib) 

= resIV {{^a - TZaY (O^ - TZb)) + resTV {Y'^Y' [Hb, Ha]) 
poly =0{x-'^) 

= TesTT(Jifjh-K^b-Knb+ KT^b +Y-^Y'[nb,na]^ 

= resTr(^ - TZ'^Qb - ^'JZb + y-^y'[Hfe,Hj) = (3.11) 

Using integration by parts (and cyclicity of the trace) on the first term here above, we obtain 
precisely the last quantity in (3.10). The Theorem is proved. Q.E.D. 
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3.1 Application to our problem 

We now apply the general definition above to our setting, with the identifications ^ = "^n, 
y = Fjv (as a formal power series at oo) and H = ^'q. We will write Y/v instead of F^v in the 
expressions below to emphasize that we consider its asymptotic expansion at oo This reduces 
the definition of the tau function to the one below 

Definition 3.2. The r-function is defined by the following PDF 

d(logrjv) = Res Tr {Y^^Y^di^o)"^^^) (3.12) 

where Yat is the formal asymptotic expansion ofV^ at infinity 





Yn = Yn\ x-'^^-'Idr^ (3.13) 

V X-^^Idd,-rj,J 

Remark 3.1. The matrix S of the previous section in our case becomes: 

/NO ^ \ 

5= {-mN-l)Idr^ (3.14) 

\0 -ruN Idd2-rM/ 

The partial derivatives of Inr^r split into two sets which have different form: 

du^ logTN = - Res TV fy^iy^/^Eii) (3.15) 

I— >oo y K J 

d,,\ogTN = Res Tr(y^iyj(^a,,(^o)^o') (3-16) 

where in the last equation the term 5t,^(^'o)^'o ^ non-zero entries only in the anti-principal 
minor of size d2 ■ 

One can notice that the situation we are looking at is a generalization of what happen in 
the one-matrix case. In the 1-matrix model, the matrix S is zero and therefore Y^ are (formal) 
Laurent series. The matrix ^'o matrix is absent in that case since there is only one potential and 
thus one recovers the usual definition of isomonodromic tau function (see [12]). Note also that in 
the derivation with respect to vj we have obtained the second equality using the block diagonal 
structure of ^'o (first row/column does not play a role). It is remarkable that the two systems 
are completely decoupled, i.e. that in the first one the matrix ^'o (containing all the dependance 
in V2) disappears and that in the second one the matrix ^0 (containing the potential Vi) also 
disappears. 
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3.2 Discrete Schlesinger transformation: Tau-function quotient 

In this section we investigate the relationship between the tau-function of Def. 3.2 and the 
partition function of the matrix modeL 

We anticipate that the two object turn out to be the same (up to a nonzero factor that wih 
be exphcitly computed, Thm. 3.4): the proof rehes on two steps, the first of which we prepare 
in this section. These are 

• proving that they satisy the same recurrence relation 

• identifying the initial conditions for the recurrence relation. 

We start by investigating the relationship between tn and tn+i', this analysis is essentially 
identical to the theory developed in [21] and used in [6], but we report it here for the convenience 
of the reader. 

From the fact that the ^'at has constant jumps, we deduce that ^'jv+i^^"'^ is an entire 
function. Moreover asymptotically it looks like: 




2— f'iV+l 




I X 











\ 





















x^^ 







\o 










-rN/ 



Yn (3.18) 



Thus, remembering that Ijv is a series x ^, Liouville's theorem states that ^'iv+i^'^ is a 
polynomial of degree one, and hence, for some constant matrices i2^,-R]y we must have 

^n+i^n^ = Rn{x) = R% + xR]^ (3.19) 

From the fact that det{Ri\f) does not depend on x (last property Thm. 2.3), we know that 
i?^^(x) is a polynomial of degree at most one as well (this is easy if one consider the expression 
of the inverse of a matrix using the co-matrix) . 

Comparing the asymptotics of ^'at+i and i?Ar(x)^Ar term-by-term in the expansion in inverse 
powers of x and after some elementary algebra one obtains ([20] Appendix A): 

Rn{x) = EagX + Rn,o and Rjj^{x) = Eix + R~^q (3.20) 
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Here we have introduced the notation ao = r^ + l which corresponds to the index of the column 
where the coefficient is to be found in the asymptotic of *I'7v+i*I'jv^- This notation is the 
standard notation used originally by Jimbo-Miwa in a Schlesinger transformation. The matrix 
{RNfi)a,f3 is given by: 

" = "0 fvt \ , -(^7V,l)ao,l -(^N,l)ao,/3 



(^JV,l)aO'l 

a = 1 

JV,1 jao,l 

and (-R^|o)a./3 is given by: 

/3 = ao = ! /3/cKo,l 
a = ao (i^Af,i)Qo,i 
a — J- — 7v — ^ 7v — S — r(JAr,i)i,i —jv — 

(i^JV.ljtiQ,! (yN,l)ciQ,l ^ iv,j./±,i (,^]V,ljao,l 



(3.21) 



(3.22) 



While the formulae above might seem complicated, we will use the two important observa- 
tions: 

EaoRj^fl + RNflEi = Rj^^qEco + EiRna = (3.23) 

R]^/'{x)R'j^{x) = Rjf^QEao does not depend on x. 

The recurrence relation satisfied by the sequence {tn} is derived in the next theorem. 

Theorem 3.2. Up to multiplication by functions that do not depend on the isomonodromic 
parameters (i.e. independent of the potentials Vi, V2J the following identity holds 

— = (^i)i,ao (3.24) 

TN 

Proof The proof follows [21] but we report it here for convenience of the reader. Consider 
the following identity 

^iv+i = = RnYn^^o (3.25) 
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This implies that 

yjv+i = RnYn (3.26) 



Taking the derivative with respect to x gives: 
Therefore we have: 



■^AT+l^iV+l — R'n^n + Yn'^Yn (3.27) 



dlogTAT+i - dlogrjv = Bbs^Tt{{Y^^R],^R'^Yn + Y^^Yj^ - Y^^YN)d{^o)^^') 

= R^Tv{Y-^R-^R'^YNd{^o)%^) (3.28) 

X— »oo 

We now need to "transfer" the exterior derivative from *o to Y^. This can be done using that 
^ = ljv*0) so that 

N 

d* = d(yiv)*o + YNdi^o) 

N 

Equivalently: 

YATd^o^o '^w' = - dY^Y^^ (3.29) 

N N 

Inserting these identities in the tau quotient we obtain the relation 

dlogrjv+i - dlogTAT = Res TV ( Rt} R'j^d{qf)^-^ - Rl^^ R'^cIYnY-A (3.30) 

x-*oo \ N N J 

The first term is residueless at oo since d^^~^ is polynomial in x and RZ^R'^ does not 

NN 

depend on x. Therefore we are left only with: 

dlogTN+i - d log rjv = - Res Tt{R]^^ R'^dYNY^^^) (3.31) 
A direct matrix computation using the explicit form of i?jv yields 

dlogTAT+i - dlogTiv = dlog((yAr,i)i,ao) (3.32) 

and hence 

= (il)i,ao (3.33) 

The last equality is to be understood up to a multiplicative constant not depending on the 
parameters uk and vj in r. Q.E.D. 

In order to complete the first step we need to express the entry (Yi)^ in terms of the ratio 
of two consecutive partition functions. This is accomplished in the following section. 
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Theorem 3.3. For the matrix the asymptotic expansion at infinity (2.10) is such that 

(^^,l)l,ao = {vd, + ifhN = {Vd,+lf^^ (3.34) 
where S and qq G {0, 1, . . . , (i2 — 1} are defined by the following relation 

N = d2S + ao-l (3.35) 

Proof In order to compute (yAr,i)i,ao it is sufficient to compute tlie leading term of tlie 
expansion at oo appearing in tlie first row of the matrix r^r. Recalling the expression (2.12), we 
start by the following direct compuation using integration by parts 

dzdw TTN{z)z'w''-^e-^^^'^~^^^'^+''" = [[ dzdw TTNiz)e~^^^'^w''-^e"^^^'"^ ^ {e'"") 



dw'^ 

dzdu;7r^(z)<7rf,,+fc_i(u;)e-^^(^)-^^(^)+^'" (3.36) 

where qd2i+k-iiw) is a polynomial of the indicated degree whose leading coefficient is "y^^+i- '^^^ 
last RHS is if d2i + k — 1 < N because of orthogonality. If d2i + k — 1 = N the integral gives 
■J^rfj+i^^ normality conditions concerning our biorthogonal set. This computation allows 

us to expand the Cauchy transform of (rjv)i^oo ii^ar oo as follows: 

'^■mJJ^ z-w 

= // rf^<i-^(^)^---^e-^^(^)-^^(^)+^- 

i=o 

+ —^11 dzd«;^z^w;"o-ie-^^(^)-^^(^)+^"' + ©(x-^-^) 
27n x^~^^ J J ^ X — z 

(3.37) 

By orthogonality the first sum vanishes term-by-term and the leading coefficient of the second 
term is ff^+i^A^- Q.E.D. 

Recalling that the r-function is only defined up to a multiplicative constant not depending 
on N nor on the coefficients u^and Vj, we have 

[Vd2+i) — (3.38) 

where N = d2SN + o-q — 1 Hence for every uq: 

TNZn, = ^JVT„o(?^d2 + l)^^="0 ' (3.39) 
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/I 


Co(l) Ci(l) . 









1 



























... ■ 


1 





One would like to take no = because it enables explicit computations. As we will prove now 
there is a way of extending naturally all the reasoning down to 0. 

The RHP for Tj\f (Thm. 2.1) is perfectly well-defined for = and has solution 



(3.40) 



Consequently we can take 

tnZo = {vd,+i)^^=o ^'Znto (3.41) 

Also note that Zq = 1 (by definition). 

We can compute tq directly from Def. 3.2 because of the particularly simple and explicit 
expression of vJ/q = Tq^'o- 

dlnro = resTr {Y^^Yt^d'^o'i>-^) (3.42) 

We claim that this expression is identically zero (and hence we can define tq = 1); indeed, 

/O * ... *\ 
... 



v-W' 



\0 





0/ 



(3.43) 



and 



d^'o(x)^'oHa;) 





* 



\0 * ... 

so that the trace of the product is always zero (even before taking the residue). Combining the 
two results together gives the following theorem: 

Theorem 3.4. The isomodromic r-function and the partition function are related by: 



(3.44) 



N 



2^i = ^3 



where we recall that Sj is given by the decomposition of j + 1 in the Euclidian division by 

d2 



S,=E 



A short computation of the power in f^j+i gives: 
VA G N : ^jv = {Vd2+i 



TN 



where ajy = E 
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The presence of the power in Vd2+i is due to a bad normahsation of the partition function 
itself {Zj\f) and can be easily cancelled out by taking w^j+i = ^ from the start (it is just a 
normalization of the weight function). Moreover it is not surprising because in the work of 
Bergere and Eynard [2], all results concerning the partition function and its derivatives with 
respect to parameters have special cases for u^^+i and fdj+i- It also signals the fact that the 
RHP is badly defined when Vd2+i = because the contour integrals involved diverge and the 
whole setup breaks down. Indeed if w^j+i ~ this simply means that V2 is a polynomial of 
lower degree and thus the RHP that we should set up should be of smaller size from the outset. 

4 Outlook 

In this article, we have restricted ourselves to contours going from infinity to infinity. This allows 
us to use integration by parts without picking up any boundary term. A natural extension of 
this work could be to see what happens when contours end in the complex plane, and especially 
study what happens when the end points moves (models with hard edges). This generalization is 
important in the computation of the gap probabilities of the Dyson model [27], which correspond 
to a random matrix model with Gaussian potentials but with the integration restricted to 
intervals of the real axis. 
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A Factorization of 



Starting from the definition of the last d2 columns of (2.18) we observe that 

N 



2ivr ./r. 



ds 



X — z 



s — w 



'dwdz (A.l) 

= J];(f ^)„^,p(Fo)p,g(W^o)g,fc = (r^V ^0W^0)m,fe (A.2) 

This proves Thm. 2.3. 

B Bilinear concomitant as intersection number 

We recall very briefly the result of [3] stating that 



dx + dz 



-w{x)f{z] 



VM.^M^^^ir,-s)-V,M+V,is) ^ 2ivrr#f = constant . 
rJr V-s 

(B.l) 

The last identity is obtained by integration by parts and shows that the bilinear concomitant is 
just the intersection number of the (homology classes) of the contours r,r. More precisely we 
get that: 

d f f , , yji.'n) - V2{s) x(r,-s)-V2{v)+V2{s) 



dsdr] 

dx ./r Jr Tj — s 

dsdr]{yi{ri) - y2'(s))e^(''~")-^2('')+^2(^) 

ds(ir/^(-e-^2(»?))gX,,g-xs+V2(s) _ f f^^^d^^V2{s)y-xs^xr,-V2{r,) 

T drj Jy Jr ds 

= X j y^(is(irye^''-^"-^2(r,)+y2(s) _^ j dsd7/e^''-^"-^2(>7)+y2(s) 

= ^ ^ ^ (B.2) 

The matrix expression shows that the pairing is indeed a duality since the determinant is nonzero. 
The undressing matrix ^'q (that was originally introduced in Thm. 2.4) is thus 





" 1 












Vd2 + l 


*0 = 




Vd2 
Vd2 + 1 


Vd2 + l 


Vd2 + l 



1 



/i 



/l /2 
/l /2 

(^2-1) 



fd2 
fd2 



Ad2'l) 
Jd2 



(B.3) 
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where the Wronskian subblock in the second term is constructed by choosing 6,2 homologically 
independent contour classes for the integrations f ; 

Ux) := [ e-^^+^^(^)ds , k = l,...,d2. (B.4) 

The dressing matrix *o exhibits a Stokes' phenomenon (of Airy's type) which is the inevitable 
drawback of removing the .x-dcpcndcncc from the jump matrix. We can now compute the jumps 
and see that it does not depend on x. For the k-th column we have: 

(B.5) 

gives: 

V'L^)(x)+ = V'^^)(x)_ + Vk°)(x)y" j dsdw ^^^'l ~_ ^2(^) e-v.(^)+y.(.)+x(^-.) (B_6) 
= ^(,^)(x)_ + Vk°)(x) ^ x,,-(rf #f ,) , (B.7) 
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